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^ A systematic LSDA+U study of doping effects on the electronic and structural properties of single 

layer C0O2 is presented. Undoped C0O2 is a charge transfer insulator within LSDA+U and a metal 
with a high density of states (DOS) at the Fermi level within LSDA. (Co02)^"~, on the other hand, 
0^ ' is a band insulator with a gap of 2.2 eV. Systems with fractional doping are metals if no charge 

orderings are present. Due to the strong interaction between the doped electron and other correlated 
Co d electrons, the calculated electronic structure of (Co02)^~ depends sensitively on the doping 
' ' level X. Zone center optical phonon energies are calculated under the frozen phonon approximation 

Q , and are in good agreement with measured values. Softening of the Eg phonon at doping x ~ 0.25 

■ seems to indicate a strong electron-phonon coupling in this system. Possible intemediate spin states 

I ' of Co ions, Na ordering, as well as magnetic and charge orderings in this system are also discussed. 

Oh' 
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PACS numbers: 71.20.-b,71.27.+a,74.25.Kc,75.20.Hr 



C3 . I. INTRODUCTION 



The recent discovery of superconductivity in hydrated ^a.xCoO-^r' has generated renewed interest in this technolog- 
ically important material. Na2:Co02 has been known for several years as a potential thermoelectric material which 
Q > exhibits an unexpectedly large thermoelectric power and at the same time a low resistivity^. Although the origin of 
' the large thermopower remains a subject of active investigation, strong correlations between Co d electrons and spin 
entropy are believed to play a critical rol o^i'^i^ . Beside its unusual thermoelectric properties, Naa;Co02 {x in the range 
0.5 ~ 0.75) is also known for having a Curie- Weiss type of susceptibility instead of a Pauli paramagnetic behavioriSi^, 
^ ^ which would be more compatible with its metallic conduction^. Although there have been reports of a weak magnetic 
ordering transition in Nao.75Co02 at r,„ = 22 KSiS, no such transition has been observed down to 2 K for systems 
with lower Na contents. Compared to the vast experimental work that has been done on this material, theoretical 
study seems to have lagged behind. 
^ First principle calculations of electronic and magnetic properties of strongly correlated systems such as Naa;Co02 
~l have always been a challenge. Although the local spin density approximation (LSDA) to the density functional theory 
(DPT) has been applied to various systems with great success, it is well known that the LSDA fails in many aspects 
"p^i , when applied to late transition metal oxides in which strong correlations between d electrons play an important role. 

■ Por example, LSDA fails to reproduce the insulating, antiferromagnetic (APM) ground state for several transition 
^ \ metal oxides, including the parent materials of high transition temperature (Tc) superconductors. Not surprisingly, 
I ■ there have been several attempts to improve the L(S)DA to take into account (at least partially) the strong electron- 
' O ' electron interactions in first-principle calculations. One of the simplest, yet very successful schemes, was proposed 
^ I by Anisimov et aljiS^ii^iSii^: the LDA+U method. In this paper, we present a systematic LSDA+U study of doping 
effects on the electronic and structural properties of Naa;Co02 (0 < a; < 1) using a recently implemented rotationally 
invariant LSDA+U methodic within the pseudopotential plane-wave formalist 



sJ . II. THE LSDA+U METHOD 

a, 

It is now well established that the failure of LSDA, when applied to late transition metal oxides, is largely due to 
an insufficient treatment of on-site Coulomb interactions between the rather localized d electrons. LSDA attempts to 
account for the Coulomb interaction via an averaged potential depending only on local spin densities. Consequently, 
magnetic moment formation is driven mainly by the spin polarization energy within LSDA. Orbital polarizations, on 
the other hand, play less important roles in LSDA and the occupation of localized orbitals does not depend sensitively 
on their orientation (symmetry). Por highly localized electrons, however, the Coulomb interaction should be applied 
to the localized orbital as a whole and is better described by the Hubbard or Hartrec-Pock (HP) type of theory. The 
formation of a local moment is therefore a result of both spin and orbital polarizations. Unfortunately, although 
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the HF aproximation and other correlated quantum chemistry methods have been apphed to atomic and molecular 
systems with tremendous success, their application to solids has been limited. Screening effects, which are relatively 
weak in atoms and molecules, are usually important in solids and cannot be easily included within the HF theory. 
Higher level quantum chemistry calculations, on the other hand, are computationally too expensive for most solids. 

The LSDA+U method attempts to incorporate the orbital specific screened Coulomb interaction while retaining 
the simplicity of LSDA. In LSDA+U, the energy fmictional consists of three contributionsifiiii: 

^LSDA+U[^<T(j.)j^^|] ^ E^^^'^[p''{r)]+E^[{n''}]- E'^'^iin"}], (1) 

where i<;LSDA ^j^g usual LSDA energy functional for spin densities p"^ (cr =tii); is a Hubbard or HF type 
of interaction arising from localized electrons (described by orbital occupation density matrices n"^) and E'^'^ is a 
"double-counting" term to be defined later. In the rotationally invariant LSDA+U methodic, E'^ takes the familiar 
HF form: 

E" = \ Y1 (('^i''^2|^^''1m3,TO4) - W(mi,m2|U'=''|m4,m3))<„^ ,„^, (2) 

{m},{a} 

where the matrix elements of the screened electron-electron interaction V'^'^ can be expressed approximately as a sum 
of Slater integrals F'': 



(mi,m2|U'"'|TO3,TO4) 
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jRUri)RKr2)dridr2, 



(4) 



An ^ 

ak{mi, 1713, 1712, 1114) ^ 2fc+ 1 E O^lmil^kqlYlms) ■ (5) 

q— — k 

Here e is the dielectric constant of the system, r< and r> are the smallest and largest values of ri and r2, and Ri 
is the radial wavefunction of the localized electron. It should be pointed out that the Slater integrals are not well 
defined in solids and expression ^ is only an approximation. For d electrons (/ = 2), three Slater integrals, 
and F*, are needed. The Slater integrals relate to the familiar Coulomb (or Hubbard) U and exchange J parameters 
a, as U — F^ and J — {F^ + i^'*)/14. A further simplication can be achieved by the observation that F^ / F'^ « 0.625 
for most c?-electron systemsii^i^. The Slater integrals F^''s (or equivalently, U and J) are fixed parameters in our 
calculations. In principle, however, one could calculate these parameters self-consistently. The double counting term 

j^,, n{n — 1) ,71^ (n^ — 1) n-'-(n-^ — 1), , , 

E^^[{n^}] = U \ > - J[^^ + (6) 

is the averaged electron-electron interaction already included in LSDA, assuming that LSDA gives the overall Coulomb 
and exchange energies correctly. The double counting term is not uniquely defined and there have been some discus- 
sions in the literature concerning other possible forms and their effects on the calculated properties^. In the above 
expressions, 71'^ =Tr(n'^) and n = + , while the density matrix n'^ for localized (e.g., d or /) orbitals remains 
to be defined. Identifying localized orbitals is trivial in computational methods using atomic basis sets such as the 
linear muffin-tin orbital method (LMTO). In the pseusopotential plane- wave method, this is less obvious and can be 
done by projecting the wavefunctions 5"^^. onto pseudoatomic orbitals Ri{r) calculated with an appropriate atomic 
configuration (3d^4s° for Co in this work): 

<«i,™. =E(^"il*»fc)(*"fcl"^2). (7) 

n,k 

We have used the abbreviation |m) = \RiYi„i). Note that proper symmetrizations of the density matrix are needed if 
the above summation is carried out in the irreducible Brillouin zone (BZ). Diagonalization of the density matrix then 
gives the symmetry-adapted local orbitals and their occupation numbers. 
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Applying the variational principle to the energy functionals defined in Eqs. (0, lO and lO, we have, in addition 
to the usual one-electron LSDA Hamiltonian, an orbital-dependent correction term 

SV^= ^ |mi)^C,,™.("i2|, (8) 

mi ,1712 

where the matrix elements 

7113 ,7n4,cr' 

The resulting one-electron problem 

{HlsBA + s^niKk) = KklKk) (10) 

can be solved self-consistently. Due to the presence of the orbital-dependent potential V , it is more convenient to 
solve Eq. (|10|l in two stepsi^. First, we solve an auxiliary LSDA problem 

HlsBAlKk) = <k\Kk) (11) 

to obtain an orthogonal basis {<&J^j,} and the corresponding eigenvalues {e^j,}- Note that the one-electron Hamiltonian 
^LSDA constructed using the electron density determined by minimizing the LSDA+U (not the LSDA) energy 
functional. We then construct the full LSDA-I-U Hamiltonian matrix 

[^fLSDA+u(fc)]n,n' = Sn,n<k + {KklSVl^n'k) (12) 

on the subspace of interest (e.g., Co d orbitals) and diagonalize it to obtain the LSDA-fU wavefunctions 5'^^, and 
band energies The occupation matrix n*^ and charge density p'^ are then constructed for the next iteration until 
the self-consistency is achieved, while fixing the parameters U and J. 



III. COMPUTATIONAL DETAILS 



Na2,Co02 assumes a layered structure in which C0O2 and Na layers alternate along the c axis. The electronically 
active C0O2 layer consists of edge sharing CoOe octahedra with magnetic (Co) ions forming a frustrated triangular 
lattice (see Fig. ^ . The oxygen octahedra are distorted considerably - compressed along the body-diagonal direction 
of the embedding rocksalt structure and stretched in the perpendicular plane. The distortion presumably depends on 
the doping level (as will be shown later). This active C0O2 layer is believed to be responsible for various abnormal 
electronic, magnetic and transport properties of the system and is the focus of the present study. The sodium layer 
is disordered, with Na ions distributed among two distinct, partially occupied sites (Wyckoff indices 2b and 2d). In 
order to avoid inconvenient (i.e., large) unit cells for systems with fractional doping, the effects of Na is modeled 
by corresponding electron doping into the C0O2 layer in our calculations. The excess electrons are then balanced 
by a uniform postive background. In real systems, the presence of Na potentials and small strains associated with 
them, as well as interlayer interactions, might have some additional effects on the electronic properties. We believe, 
however, our model captures the essential physics and the effects of Na ions are minor if not negligible, as will be 
discussed later. (Alternatively, one may employ the virtual crystal technique to overcome the large unit cell problem.) 
Our treatment might become even more exact in the case of hydrated compounds since water molecules are likely to 
screen out the Na potentials. Of cource, the interaction between H2O and C0O2 layers is another issue that deserves 
further investigations. In our calculations for single layer C0O2, we fix the in-plane lattice constants a = b = 2.823A-ii, 
regardless of the doping level. Small variation in lattice constants should have negligible effects on our results. The 
separation bewteen layers is set at 6.5 A to ensure no significant interlayer interactions. Therefore, the only structural 
parameter allowed to relax is the oxygen z coordinate, which turns out to be rather sensitive to the doping level, as 
will be discussed later. 

The above assumptions significantly simplify our calculations. There are, however, other difficulties arised in 
studying correlated systems due to the existence of several competing charge and/or magnetic orderings. Since the 
magnetic and/or charge ordering energies are usually very small, it is sometimes difficult to distinguish between 
different ordering states based only on their energy difference. This is particularly true in magnetically frustrated 
systems such as C0O2, as is evidenced by experiments where no long range magnetic ordering is observed for Naa:Co02 
except for x ^ 0.75^^. Nevertheless, it was suggested that a short range ferromagnetic (FM) ordering might be 
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preferred in these systemsi^. Therefore, we will primary concentrate on the ferromagnetic (paramegtic for x = 1) 
phase in this paper. Different orderings will be discussed briefly. 

We employ the LSDA+U method as described in section II to study doping effects on the electronic, structural 
and magnetic properties of single layer (Co02)^~ (0 < x < 1). The fc-point set is generated by the Monkhorst-Pack 
scheme^^ with a density of 12x 12x2. The plane-wave energy cutoff is set at 250 Ry to ensure the convergence of the 
calculations. Such a high plane- wave energy cutoff is necessary for systems containing very localized d electrons using 
norm-conserving pseudopotentials^fliSi. Since there has been no theoretical and/or experimental determination of U 
and J for Co d electrons in Naa;Co02, we adopt a moderate U = 5.5 eV and a, J — 0.9 eV in our calculations and 
neglect their doping dependence. The exchange parameter J is of the order of 1 eV for most later transition metal 
oxidesiS and Singh has given an estimate oi U = 5 ^ 8 eV for Naj^CoO^i^. Similar values of U (5.4 and 5.0 eV) have 
been used in previous studies on this systemSSiS^. In general, LSDA-I-U results are insensitive to small variation of 
these parameters. 

IV. RESULTS 

A. Electronic structure of C0O2 and CoOj °~ 

We first study the undoped parents material C0O2 in its ferromagnetic phase. Single layer C0O2 has a D^^i point 
group symmetry, which derives from the cubic {Oh) symmetry after a distortion along the [111] direction. The Co 3d 
orbitals split into a triplet (^29) and a doublet (cg) under the influence of the octahedral (cubic) crystal field. Upon 
further lowering the symmetry, the t2g states split into Cg and aig levels. The t2g derived eg states then mix with the 

original Cg ones, forming two new doublets e^^-* and el^\ Of course, the degree of this mixing increases with increasing 
trigonal distortion and Co d derived eg states will further hybridize with the O p states. It is generally believed that 
the relevant low-energy electronic states of C0O2 are predominately of Co d character and can be interpreted in terms 
of those of the Co ion. For example, the electronic structure of undoped C0O2 in its low spin state (S* = i) may be 

understood in terms of Co''"'"(eJajge^). Thus upon electron doping, it becomes a doped spin-i system. However, due 
to the strong mixing between O p and Co d states in these systems, the validity of such a simplified picture needs to 
be carefully examined. 

Figure[21compares the LSDA and LSDA+U density of states (DOS) of Co02. The undoped parent material C0O2 is 
a Mott-Hubbard insulator (or charge transfer insulator according to the ZAS classification^'*) as predicted by LSDA-I-U. 
In contrast, our LSDA calculation gives a metallic ground state with a rather high DOS (2.5 electrons/eV/cell) at 
the Fermi level, which is consistent with the LSDA result of Singh^* . The local spin moment of Co calculated within 
LSDA-I-U is about 1 /is, as expected for a spin-i system. LSDA, on the other hand, gives a local moment of 0.58 hb 
due to lack of orbital polarizations. The octahedral crystal field splitting (^ 3.0 eV) of Co d orbitals, i.e., the splitting 
between the occupied t2g and the unoccupied eg states of the d-orbitals, can easily be estimated from the LSDA 
results. Further splitting of the t2g states (the separation between the two peaks within the triplet labelled t2g{d + 6p) 
in the upper panel of Fig. O due to the trigonal distortion is about 1.0 eV within LSDA. Determining these values 
from LSDA-fU results, however, is more involved since additional splitting due to the Coulomb U can not be easily 
decoupled. The explicit removal of the self-interaction in the screened HF interaction term in the LSDA-I-U method 
separates the occupied and the unoccupied d states and pushes the occupied ones below the O p levels. Consequently, 
the top of the valene band has predominately O p character within LSDA-I-U, contrary to the Co d character in LSDA. 
Another interesting observation is that the hybridization between the occupied Co d and O p states is enhanced within 
the LSDA-I-U method. Whereas the DOS calculated with LSDA clearly shows p-dominate t2g and eg states (labelled 
t2g{p + Sd) and eg{p + Sd) in Fig. [^J and d-dominate t2g states (labelled t2g{d + Sp)), this separation between d and 
p states becomes less obvious in LSDA-I-U results. Therefore, our results suggest that both Co d and O p states need 
to be considered when one attempts to derive an effective low energy model hamiltonian for this system. 

To better understand the electronic structure of C0O2, we further project the wavefunctions onto the symmetry- 
adapted Co d orbitals (i.e., the eigenfunctions of the density matrix defined in Eq. {T))), as shown in Fig. |31 The 
corresponding orbital occupations are given in Table As we have mentioend above, the t2g triplet splits into eg 
and aig under the influence of triangular crystal field. Although this splitting (~ 1.0 eV) is insignificant in LSDA, 
the strong on-site Coulomb interaction included in the LSDA-I-U pushes the minority spin aig up so that it becomes 
completely unoccupied. The separation between the occupied and the unoccupied aig states is [/ -I- J ~ 6.4 eV, as 
shown in the upper panel of Fig.|21 The t2g derived eg doublet, labelled ef\ is nearly fully occupied as expected (see 

Table nj. For the other doublet, i.e., eg^\ one would expect it to be unoccupied from the simple molecular orbital 
analysis. However, this does not seem to be the case in our calculations due to strong mixing between Co d and O p 
states with compatible symmetry. We call these hybridized doublet eg(pd). Both the occupied and the unoccupied 
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d orbital occupation 







LSDA 


LSDA+U 






0.845 


0.926 


Majority 


^9 


0.907 


0.919 


Spin 


^9 


0.457 


0.506 




Total 


3.57 


3.78 




aig 


0.668 


0.097 


Minority 




0.725 


0.916 


Spin 




0.434 


0.383 




Total 


2.99 


2.69 


Total d electrons 


6.56 


6.47 


Co Spin Moment (ps) 


0.58 


1.09 



TABLE I: Comparison of d electron occupation of undoped C0O2 between LSDA and LSDA+U results. For doublets, the 
number of d electrons is the occupation times two. 



eg{pd) have nearly equal O p and Co d characters. As a result, the valency of the Co ion in C0O2 deviates substantially 
from its norminal value 4+. This indicates a coexistence of ionic and covalent bondings in this system: Whereas the 
Co 4s electrons are fully ionized, the d electrons are better described as covalently bonded with O p electrons due to 
the significant overlap between the two wavefunctions. 

Upon doping one electron to the C0O2 layer, the system becomes a non-magnetic insulator with a band gap of 
about 2.2 eV (see Fig.0J, which compares favourably with the measured band gap (~ 2.7 eV) of a similar material 
LiCoOz^. Interestingly, LSDA also predicts a non-magnetic insulating ground state but with a smaller gap of about 
0.8 eV (see Fig. 0J. This is expected since the crystal field splitting of Co d states within LSDA is larger than the 
bandwidth of the t2g and Cg sub-bands. Therefore, NaCo02 is a band insulator within LDA. The on-site Coulomb 
interaction thus contributes about 1.4 eV to the calculated band gap, which is only a fraction of the parameter U . This 
is an indication that the results of LSDA-I-U calculations are not especially sensitive to the parameters. Although both 
LSDA and LSDA-I-U give qualitatively the same nonmagnetic insulating ground state, there are significant differences 
in the calculated DOS with the two methods, especially for the occupied states. As for the undoped case, LSDA-I-U 
enhances the hybridization between the Co d and O p states. The top valence band triplet has predominately Co d 
character in the absence of local Coulomb interactions but strongly hybridizes with O p states within LSDA-I-U. Also, 
there is a small but not negligible gap between the d-dominate and p-dominate states within LSDA, a feature that does 
not exist within LSDA-I-U and was not observed in photo-emission experiments^^. Overall, our results show that the 
(i-states mostly concentrate within a 2 eV window below the valence band maximum and the oxygen p states spread 
from —7.0 to —2.0 eV. This is in good agreement with the resonant photoemission experiment results for LiCoOa^ 
where the sharp peak around -1.4 eV is assigned to the Co d final states whereas the broad structure at -5±2 eV is 
attributed to O p states. Of course, one must be cautious when comparing the calculated DOS with the photoemission 
results since both matrix element effects and correlations may change the lineshape of the photoemission spectra. 

The d orbital occupation of both majority and minority spins are significantly affected by the electron doping (see 
Tables III and HHl . Whereas the minority spin occupation increases by 0.68 electrons within LSDA-I-U, the majority 
spin occupation actually decreases by 0.41. This is not a surprising result for a strongly correlated system: Although 
the lower and upper Hubbard bands (LHB and UHB) are energetically separated, they are intimately correlated and 
are both affected by doping. Strong doping dependence of both the band energies and the spectra weights of the 
correlated d bands can also be seen by comparing the projected DOS between the undoped and doped cases (see the 
lower panel of Fig. |2|and the upper panel of Fig. 01 . Due to their mixing with the correlated d states, O p states 
are also affected by the doping. Overall, the total d electron occupation only increases by 0.27 electrons upon doping 
one electron to the system. This suggests that a substantial portion of the doped charge actually go to oxygen sites 
and that the Co valency is rather insensitive to the doping level in this particular system. This insensitivity is also 
recently reported in Na:rCoi_yMny02 systems, where Co ions are partially substituted by Mn2&. 

To ensure the validity of replacing Na ions with a uniform positive background, which reduces the complex system 
to a single-layer one, we have also done a calculation for a realistic system, i.e., a double-layer NaCo02. Fig. |5| 
compares the DOS of the model system (0002)^"" with that of NaCo02. The negligible differences between the 
DOS of the two systems suggest that our model system is appropriate and interlayer coupling is not very important in 
this system, assuming no charge and/or magnetic orderings are involved. The possibility of a subtle interplay between 
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d orbital occupation/spin 





LSDA 


LSDA+U 


aig 


0.887 


0.908 




0.889 


0.908 


^9 


0.363 


0.321 


Total 


3.39 


3.37 



TABLE II: d electron occupation of (Co02)^ '' calculated with LSDA and LSDA+U. The system is non-magnetic. For doublets, 
the number of d electrons is the occupation times two. 



Na potentials and charge orderings in the C0O2 layer will be discussed later. 

B. Electronic structure of (Co02)''~ (0 < x < 1) 

Having discussed the electronic properties of the two extreme cases, it would be interesting to see how the electronic 
structure of (Co02)^~ evolves as the doping level x varies. Fig. Elshows the evolution of the DOS as the doping level 
X changes from 0.0 to 1.0. As we have discussed, the undoped parent material C0O2 is a charge transfer insulator 
with a gap between the O p and Co d (aig) states. Upon electron doping, the originally unoccupied aig state, which 
is split-off from the t2g triplet due to on-site Coulomb interactions, becomes partially occupied and moves toward 
lower energy, touching the O p-dominate valence band at a; '--^ 0.3 and eventually merging into the rest of the valence 
band. At a; = 1.0 the t2g triplet is recovered and occupied by both spin up and spin down electrons. At the same 
time, a new gap (larger than that of the undoped case) appears between the t2g and eg states. All other occupied d 
states, in contrast, are pushed upward as doping level increases. As a result, significant spectra weight moves across 
the entire valence bandwidth due to the correlation between the doped electron and other occupied d states. 

The lowering of the energy of the minority spin aig (relative to other d states) upon electon doping is intriguing 
and should be studied in more detail. Table imi shows the d electron occupations for all symmetrized d orbitals. The 
occupation of all other d states except the minority spin aig decreases as doping level increases, with the majority 

spin eg^-* state being affected the most due to less hybridization between this state and the O p states as doping 
increases. Consequently, the effective Coulomb repulsion "felt" by the minority spin aig electron is actually reduced 
with increasing doping level. Note also that the occupation of this minority spin aig state is approximately the same 
as the doping level. However, this does not mean that all doped electrons go to Co site. In fact, it seems that 
a significant portion of the doped charge actually goes to oxygen sites, as we discussed previously for the case of 
(Co02)^ °^. Table ITTll also gives the local spin moment on Co site as a function of doping level, which decreases 
monotonically with increasing doping (m^ ~ 1.0 — a; fJ-s)- We should point out that the occupation analysis in our 
calculations is only approximate and different charge analyses may give slightly different results. The subtle changes 
to the Cg and majority spin aig occupation might be partially due to orbital and/or structural relaxation effects. 

Another factor that contributes to the doping dependence of the energy of the minority spin aig state comes from 
the doping dependence of the trigonal distortion: The crystal field splitting within the t2g triplet decreases with 
decreasing distortion. As we will show in the following, trigonal distortion increases as doping level decreases, and so 
does the splitting. 

The minority spin aig state deserves particular attention since, upon electron doping, it determines the low-energy 
electronic properties of the system and has been the subject of intensive discussion concerning its connection with the 
superconductivity observed in Nao.3Co02-yH2032i22i224^2i2i. Fig.[7|shows the band structure of (Co02)°'^^~ with the 
minority spin aig band highlighted. Interestingly, apart from a constant shift, this band can be fairly well fitted by 
a simple tight-binding model with a nearest hopping parameter t = —0.155 eV. Note that here t is not the hopping 
element between Co and O sites but the effective hopping between neighboring Co sites. The total bandwidth is 
thus 9\t\ = 1.4 eV. Since we use a on-site Coulomb interaction U = 5.5 eV, the effective superexchange between two 
neighboring Co ions is J = At^ /U = 16 meV, in reasonable agreement with a previous estimate^. Due to the particular 
dispersion and doping dependence of this band, there is a strong doping dependence to the Fermi-surface properties 
and the DOS at the Fermi level, as shown in Fig.il The DOS(£;/) increases sharply with increasing doping level for 
X < 0.1, reaches a maximum at about x ~ 0.2, and then decreases with increasing doping. The narrow doping range 
0.1 < X < 0.3 beyond which DOS(i?/) decreases rapidly with increasing or decreasing doping level is closely related 
to the Fermi surface structure of the system: At very low dopings, small electron pockets appear around the corners 
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doping 


0.0 


0.25 


0.5 


0.75 


1.0 




aig 


0.926 


0.922 


0.917 


0.912 


0.908 


Majority 




0.919 


0.916 


0.914 


0.911 


0.908 


Spin 




0.506 


0.436 


0.381 


0.345 


0.321 




total 


3.78 


3.63 


3.51 


3.42 


3.37 






0.097 


0.341 


0.566 


0.754 


0.908 


Minority 




0.916 


0.913 


0.910 


0.908 


0.908 


Spin 


"-s 


0.383 


0.379 


0.367 


0.346 


0.321 




total 


2.69 


2.93 


3.12 


3.26 


3.37 


Co Spin Moment {/j-b) 


1.08 


0.70 


0.39 


0.16 


0.0 



TABLE III: d electron occupation and local spin moment of (Co02)^ as a function of doping level. For doublets, the number 
of d electrons is the occupation times two. 



of the BZ. (Note that this metallic state may not be stable against charge orderings at very low doping levels.) The 
Fermi surface quickly extends and then shrinks with increasing doping. Compared with LSD A results where a large 
Fermi surface, as well as small pockets of holes, are predicted, there is only one large Fermi surface in our calculation, 
which agrees well with a recent experiment^. Although it would be interesting to connect this observation with 
the fact that superconductivity occurs only in a very narrow doping range (x ~ 0.3 in Naa;Co02-yH2Oi*^), further 
investigation on this subject is required. If the superconductivity in this system is of phonon origin, then the high 
DOS at the Fermi level is definitely an important factor in determining the superconducting transition temperature. 

C. Doping effects on the structural properties and possible spin-phonon interactions 

Structural properties are in general not particularly sensitive to the doping level. As we have mentioned above, we 
fix the lattice constants in our calculations but allow the oxygen atoms to relax. Fig. IHl shows the doping-dependent 
oxygen z coordinate as measured from the Co plane. The calculated O z coordinate at x ~ 0.30 (1.72 a.u.) falls 
within the measured values (1.67 ~ 1.77 a.u.) for Naa;Co02 • ?/H20ii24 but is smaller than those (~ 1.83 a.u.) for 
unhydrated systemsSl. This is reasonable since the single layer system in our calculations should mimic the hydrated 
system better than the unhydrated one. The Na potential in the unhydrated system is likely to attract the negatively 
charged O ions away from the Co layer, making the O z coordinate larger. Overall, the distance between the oxygen 
and cobalt layers expands quadratically with increasing doping level. 

Using the LSDA+U total energy functional defined in Eq. we can calculate phonon energies under the frozen 
phonon approximation^^. Single layer C0O2 (assuming paramagnetic/ferromagnetic ordering) has four zone-center 
optical phonon modes. Two of them relate to the in-plane and out-of-plane motion of oxygen atoms {Eg and Aig), 
which are Raman active. The other two (i?„ and ^2m) involve cobalt moving against oxygen and are infrared (IR) 
active. Table Hvl lists the calculated zone center phonon energies for doping levels x = 0.0, 0.25, 0.5, 0.75 and 1.0. 
The calculated energies are in good agreement with available measurements^^iSSiSi. For example, the measured Aig 
phonon energy (ranging from 71.2 to 74.4 meV^S*^ depending on the doping level and sample conditions) agrees 
well with the calculated ones (from 72.6 to 73.0 meV for doping level 0.25 < x < 0.75). The theoretical Eg phonon 
energies are 63.2 meV for x = 0.25 and 65.6 meV for x — 0.75, to be compared with the measured values 56.8 ~ 61.2 
meV— . The measured energy for the IR active mode i?„ is about 70.7 meV for Nao.srCoOz^, which also compares 
favorably with our theoretical value (74.7 meV for x — 0.5). There has been no measurement for the mode 
so far. In general, the zone center optical phonons are not sensitive to the doping level. There is, however, one 
interesting exception: At ~ 0.25 doping, the Eg phonon softens, decreasing from 67.8 meV for x = 0.0 to 63.2 meV 
for X = 0.25, and has significant anharmonicity. (The calculated harmonic phonon energy is only 56.5 meV.) The 
softening of this phonon mode at doping level x < 0.2 seems to indicate a strong electron-phonon coupling and we 
believe that the high DOS at the Fermi level for x ~ 0.25, together with strong electron-phonon couplings between 
this mode and conducting states, is responsible for the phonon softening and anharmonicity and may eventually lead 
to a superconducting phase transition. 

In magnetic systems, phonons might interact with the spin degree of freedom. In fact, we observe a strong correlation 
between the O z coordinate (relates to the Aig phonon displacement) and the local spin moment on Co sites (see Fig. 
I10|l . No such correlations were found for other phonon modes. This raises the possibility of interactions between Aig 
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zone center phonon energy (meV) 


doping level 


0.0 0.25 0.5 0.75 1.0 


Ala 


72.3 72.6 73.4 73.0 71.4 


A2u 


72.8 73.6 74.0 77.0 80.0 


Ea 


67.8 63.2 (56.5) 64.6 65.2 62.7 


Eu 


75.3 74.9 74.7 74.2 72.7 



TABLE IV: Doping effects on zone center phonon energies calculated using frozen phonon approximation. All phonon modes 
are fairly harmonic except the Eg mode in the case of doping level x ~ 0.25. For x — 0.25, the calculated harmonic frequency 
is shown in parentheses. 



phonons and magnons in this system. 

V. DISCUSSIONS 
A. Other possible spin configurations 

So far we have assumed a low spin state for Co ions and the doped electrons always go to the minority spin aig 
conduction band. Thus the local spin moment on Co decreases with increasing doping level and vanishes at x = 1.0. 
However, this does not seem to be consistent with the observations that sizable effective magnetic moments fics 
exist at doping level x ^ 0.75. For example, magnetic susceptibility measurements of Nao.75Co02 give a /JcS of 2.74 
/is /Co assuming only Co*+ ions contribute to the Curie constant CS. This is much larger than the "spin-only" 
value {g\/s{s + 1) = -s/S/is) of Co''"'" ion in its low spin state. If, on the other hand, all Co sites are assumed 
equivalent, the results yield a of 1.37 /is/Cc^, which is again not compatible with the calculated spin moment. 
This apparent discrepancy raises the possibilities of unquenched orbital moments and/or other spin configurations 
of Co ions in this system, especially for high doping levels {x > 0.5). The spin states of Co ions in C0O2 are 
determined by several competing factors such as the crystal field splitting, Hund's rule coupling and the screened 
on-site Coulomb interactions. Therefore, depending on the relative strength of these factors, some Co ions might 
adopt an intermediate-spin state. Here we explore such a possibility. 

As we can see from Fig. 12 the energy of the unoccupied majority spin eg is only slightly higher than that of 
minority spin aig in the undoped single layer C0O2. In the real system, especially in unhydrated Na3;Co02, however, 
the situation might be more complicated. The presence of a possibly ordered Na potential might enhance or induce 
charge orderings in the C0O2 layer— and the energy level of the two unoccupied states, eg and aig, might get reversed 
on some Co sites. A smaller crystal field splitting or larger U could pull down the eg state or push up the aig 
one. The doped electron could then go to the eg majority spin state, resulting in an intermediate-spin state for the 
corresponding Co ions. Therefore, we propose the following scenario for unhydrated Naa;Co02 system: At low doping 
level, all Co ions are in their low spin states. As doping level increases, Na ions tend to order themselves to minimize 
the Coulomb energy. (We will discussed Na orderings in more detail in the following.) This ordering might then 
enhance or induce charge ordering in C0O2 layers. Due to charge ordering in the C0O2 layer, Co ions then have a 
different symmetry and chemical environment, leading to different crystal field splitting and/or on-site interactions. 
Under certain circumstance, the unoccupied majority spin eg state might be lower in energy than the aig state. The 
local moment of these Co ions will then increase with increasing doping. We have calculated a ferromagnetic system 
with Co ions in their intermediate spin states for doping level x = 0.5 and found that the local spin moment is 1.52 
/Us/Co. Not surprisingly, the energy of the intermediate spin state is higher than that of the low-spin state (by ~ 0.5 
eV/Co). However, charge orderings might reduce or inverse this energy difference, resulting in an intermediate-spin 
ground state. It is also plausible that this kind of rich degeneracy of spin states, a result of nearly perfect balance 
between the crystal field and Hund's rules effects, is responsible for the unusually high thcrmopower in this system. 

B. Na ordering, magnetic and/or charge orderings 

It is usually assumed that the Na layer is disordered in Na2:Co02 and the primary effect of Na is to provide 
electrons to the C0O2 layer. In strongly correlated systems, however, charge and/or spin orderings usually happen at 
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an extremely low energy scale and seemingly insignificant interactions can sometimes result in profound changes in the 
electronic structure. Although we have shown that the Na potential has minimal effects on the calculated electronic 
structure of ferromagnetically ordered C0O2, it is still unclear the exact role Na plays in determining the properties 
of the system, especially at high doping when the Na layer is likely ordered. Whether or not Na becomes ordered 
depends on a competition between entropic and energetic factors. If energetics dominate, an ordering might occur. 
In fact, there has been increasing experimental evidence that the Na layer might be ordered at some doping levels. 
Foo et al. observed a -s/S x 2 ordering in Nao.5Co02 and discussed the possible effects of Na ordering on its electronic 
properties'^^. Shi et al. reported a 2 x 1 superstructure in Na2:Co02 for x > 0.75 and ascribed it to a possible Na 
ordering"*". 

Since Na are fully ionized, it might be possible to discuss their ordering by simple electrostatic and entropic 
arguments. (Chemical interactions between Na and C0O2 layers might also play a minor role but more difficult to 
characterize.) At low Na concentrations, entropic effects should dominate and the Na layer is likely disordered. As 
doping level increases, however, Na ions will tend to organize themselves (at least partially and locally) to minize the 
ionic repulsion, since it costs energy to place two Na ions in neighboring 26 and 2d sites (see Fig. [TT|l . This does not 
necessarily lead to a long-range ordering since there could be many (nearly) degenerate local orderings. However, if 
there exist an ordered pattern which has significantly lower energy (compared to the thermal energy) than others, a 
long range ordering might result. For x = 0.5, we indeed find a particular arrangement of Na ions which is compatible 
with the observed \/3 x 2 ordering'^^ and has as much as 0.3 eV/Na lower Coulomb energy than other configurations 
with similar unit cells (see Fig. 1111) if only in-plane Coulomb interactions are taken into account. It might be possible 
that the interplay between the ordering in the Na layer and the charge ordering in the C0O2 layer is responsible 
for the observed insulating behavior in Nao.5Co02, as discussed by Foo el aP^ . For x = 0.75, the situation is more 
complicated. We find many possible orderings with similar energy. The lowest-energy pattern is show in Fig. 1111 with 
filled black circles denoting occupied and gray circles partially occupied (50% for x = 0.75) Na sites. Interestingly, 
this pattern of low energy ordering is also consistent with the reported Na superstructure for x > 0.75^. The low 
energy ordering pattern for a; = 0.5 is rather exclusive in the sense that further addition of Na to this structure will 
result in occupation of neighboring 2b and 2c? sites thus increase the Coulomb energy sharply. 

In general, we find it very energetically unfavorable to occupy neighboring 2b and 2d sites. At high doping (e.g., 
X > 0.75), Na ions tend to occupy only one of the two distinct sites within a given domain. The size of these 
domains presumably increases with decreasing temperature. If interactions between Na and C0O2 layers are taken 
into account, the two Na sites may not be equilvalent energetically. In fact, our calculations indicate that the 2c? site 
has about 0.1 eV/Na lower in energy. This is due to the different Coulomb repulsion between the Na and Co in the 
two configurations. Therefore, 2c? sites are more likely to be occupied, provided that no immediate neighboring 2b 
sites are already occupied. Although this differentiation between the two sites should be taken into account when 
discussing Na orderings, no changes to our conclusion for a Na ordering of a; = 0.5 is expected due to the large in-plane 
ordering energy. 

Magnetic and/or charge orderings in Na2;Co02 are other subjects of great interest. Kunes et al.^^, and Motrunich 
at alm^, discussed a possible -v/S x \/3 charge ordering for doping x ^ ^. Foo et al. reported a -^/S x 2 Na ordering 
at a; = 0.5 as we have mentioned above. NMR measurements also point to possible charge orderings in Na2:Co02 
for 0.5 < a; < 0.75^2ii^. However, there is no consensus in the literature on this matter so far. The magnetic 
ordering in Na2:Co02 is even more intriguing. Very weak magnetic ordering has been observed at T 22 K only 
for doping x ~ 0.75. Curiously enough, the measured magnetizations differ by two order of magnitute between two 
experimentsSiS. We try to explore the simplest AFM ordering with a 2 x 1 unit cell in this system. The AFM phase is 
found to be slightly lower in energy (~ 10 meV/Co) than the FM phase for doping level x = 0.3. However, this small 
difference could be beyond the accuracy of our theoretical treatment. Fig. E| compares the calculated DOS of FM 
and AFM phases of (Co02)'' '^~. Although the low energy valence states are not significantly affected, the bandwidth 
of the partially occupied aig state is greatly reduced as a result of AFM ordering. This raises the possibility of further 
supression of the aig bandwidth if more complicated orderings are present, which might help to account for the 
mysterious electron mass enhancementSi in this system. Unfortunately, due to the extremely small energy differences 
between competing ordering states, fluctuations among these states result in very slow convergence in self-consistent 
calculations. 



VI. CONCLUSION 

In conclusion, we have carried out systematic studies on the electronic, magnetic and structural properties of single 
layer (Co02)^~ using a recently implemented rotationally invariant LSDA-fU method within the pseudopotential 
plane-wave formalism. Both the undoped and one integer electron doped systems are insulators within LSDA+U, 
whereas systems with fraction doping are half-metal in the absence of charge ordering and assuming a ferromagnetic 
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phase. Calculated Fermi surface and zone center phonon energies agree well with available measurements. Possible 
intermediate spin configurations of Co ions, Na orderings, and magnetic and charge orderings in this system are also 
discussed. 

Although the pairing mechanism that leads to superconditivity in this system remains a subject of intensive inves- 
tigation, high DOS at the Fermi level at low doping levels, together with strong elcctron-phonon couplings, might be 
partially responsible for the superconducting transition in the hydrated systems. The role water molecules play in 
the superconducting transition is still unknown. One possibility is that the screening effects, which greatly supress 
the interaction between Na and C0O2 layers and possible charge orderings in the C0O2 layer, lead to a more homo- 
geneous electronic system in the C0O2 layer and ultimately favor a superconducting state over competing phases. 
A better understanding of the properties of this material requires more experimental work and thorough theoretical 
investigations. The interplay between the Na ordering and the charge/magnetic ordering in the C0O2 layer deserves 
particular attention, especially in the unhydrated system. 
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FIG. 1: Orthographic (upper) and perspective (lower) views of single layer C0O2. Dark balls are oxygen and gray ones are Co. 
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FIG. 2: Spin resolved DOS of C0O2 in its ferromagnetic phase calculated with LSDA (upper) and LSDA+U (lower). Solid 
curves show the total DOS and the dotted ones are partial DOS projected on Co d orbitals. 
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FIG. 3: Local DOS of C0O2 projected onto the symmetry-adapted Co d orbitals. The total DOS are also shown in dotted 
curves. 
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FIG. 4: Compaxison between LSDA and LSDA+U DOS of (Co02)^ "~. Solid curves are total DOS and dotted ones are 
projections onto Co d orbitals. 




FIG. 5: DOS of (Co02)^ ° (solid curve) and that of NaCo02 (dash curve). The DOS is not altered significantly by replacing 
Na ions with a uniform positive background. Also, interlayer coupling does not seem to be important in this system. 
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FIG. 6: Doping effects on the DOS of (Co02)^^ (0.0 < x < 1.0). Solid curves are total DOS and the dotted ones are projections 
onto Co d orbitals. As the doping level increases, the partially occupied minority spin aig state moves across the band gap 
while pushing other correlated d states upwards. The dashed line shows the Fermi level of metallic systems. 
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FIG. 7: Band structure of (Co02)'^'^^^ • Solid and dash curves are for minority and majority spins, respectively. The minority 
spin aig band is highlighted with a thick solid curve and the dotted curve is a nearest neighbor tight-binding fitting for this 
band with a hopping parameter t = 0.155 eV. 
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FIG. 8: Density of states of (Co02)^"~ at the Fermi level as a function of doping level x. 
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FIG. 10: Correlation between the local spin moment on Co and the oxygen z coordinate of C0O2. 
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FIG. 11: Possible Na ordering pattens. Upper: Two distinct Na sites in NaxCo02. At doping x = 1.0, Na will occupy only 
one of the two sites at low temperature. Middle: Lowest energy ordering pattern for x = 0.5. Lower: A low energy ordering 
pattern for x = 0.75. Gray circles indicate partially (50%) occupied Na sites. 
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FIG. 12: Magnetic ordering effects on the electronic structure of (Co02)°'^~. The sohd curve is for FM phase and the dotted 

curve for AFM phase. Although no significant changes to low-lying valence states are observed, the width of the aig conduction 
band is renormalized (narrowed) appreciately due to the ordering. 



